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The influence of ultra-fast laser pulses on electron transfer in molecular wires studied 

by a non-Markovian density matrix approach 
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New features of molecular wires can be observed when they are irradiated by laser fields. These 
effects can be achieved by periodically oscillating fields but also by short laser pulses. The theoretical 
foundation used for these investigations is a density matrix formalism where the full system is 
partitioned into a relevant part and a thermal fermionic bath. The derivation of a quantum master 
equation, either based on a time-convolutionless or time-convolution projection-operator approach, 
incorporates the interaction with time-dependent laser fields non-perturbatively and is valid at low 
temperatures for weak system-bath coupling. From the population dynamics the electrical current 
through the molecular wire is determined. This theory including further extensions is used for the 
determination of electron transport through molecular wires. As examples, we show computations 
of coherent destruction of tunneling in asymmetric periodically driven quantum systems, alternating 
currents and the suppression of the directed current by using a short laser pulse. 



I. INTRODUCTION 

Intra- and intermolecular electron transfer has been 
studied for several decades^*^ including its coherent 
control^. Lately the closely related field of molecu- 
lar electronics has attracted much interest, especially 
the transport through molecular wires^. Experimen- 
tal progres o 7 i 8 i 9 made in this field also spurred a large 
theoretical effort 5 .* 6 *. Many of those theoretical studies 
are based on a tight-binding model for the wire. The 
current-voltage characteristics is calculated either using 
a scattering approach- or an electron transfer scheme** 
while both formalisms have the same rootsi^. The well- 
known Landauer-Biittiker scattering formalismi&A4 pro- 
vides a method to compute the steady-state currents in 
systems on the nanoscale connecting two or more elec- 
trodes and later extensions of the formalism also deal 
with oscillating fieldsi£*i&. Many theories are based on 
the assumption of a weak coupling of the wire to the 
leads and employ a perturbation theory in this small 
paramete»i2»i2ii2i2Jl. In this aspect these approaches are 
very similar to Redfield-like approaches in the field of dis- 
sipative quantum dynamics. Furthermore, in some the- 
ories the wire is coupled to a dissipative environment to 
mimic relaxation and decoherence in the wire2ii22i2i2i 
or to determine current-induced light emission 25 . Ad- 
ditional effects are observed when the molecular wire is 
irradiated by a periodic laser fieldi 7 *^ 4 * 2 ^^ 7 .. Only very 
little is known for non-periodic driving fields 18,28 . One 
of these studies 2 ^ uses an extension of density functional 
theory to describe time-dependent transport phenomena. 

In this contribution we are also interested in a descrip- 
tion of the current through a molecular wire which is 
based on a quantum master equation. The aim is to 
extend those methods and results obtained for periodic 
laser field o 17 ! 241 26,27 to the area of non-periodic ultra-fast 
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laser pulses. To do so we borrow techniques developed 
in the area of dissipative quantum dynamics. Instead 
of coupling the system, i.e. the wire, to a bosonic heat 
bath it is coupled to a fermionic particle reservoir with 
which it can exchange particles. Externally applied op- 
tical fields can influence the dynamics in a direct man- 
ner by changing the wire part but also in an indirect 
way by influencing the wire-lead coupling. In the area of 
dissipative quantum dynamics these processes have been 
mainly studied for monochromatic laser fields 29,30,31,32,33 
but also for short laser pulses 3 ^ 5 .. While for monochro- 
matic fields one often employs Floquet states, as for the 
transport in periodically driven wiresi 7 ^^^^* 2 ^, the situa- 
tion is more complex for non-periodic laser pulses. Meier 
and Tannor— proposed a method how to use a special 
parametrization of the so-called spectral density of the 
reservoir leading to a set of coupled equations for a pri- 
mary and several auxiliary density matrices. Within this 
approach one automatically accounts for the influence of 
the external field on the dissipation operator. The tech- 
nique of Meier and Tannor is based on the time-nonlocal 
(TNL) Nakajima-Zwanzig identitji 3 ^ 7 -, while a similar 
scheme with auxiliary operators has been developed also 
for time-local (TL) quantum master equations 3 ^*— A 
similar technique, i.e. a direct decomposition of the cor- 
relation function, was introduced by Korolkov and Para- 
manov in earlier work—. 



In the current work the technique of auxiliary density 
matrices is used in the framework of molecular wires. 
The model Hamiltonian employed for this purpose is in- 
troduced in the following section. Section III shows the 
evolution of the density matrix using a decomposition of 
the spectral density while Section IV is dedicated to the 
derivation of the current equation used for the numerical 
calculations. In Section V, we will make a few comments 
on TNL versus TL theory. The formalism is tested for 
periodic and non-periodic laser fields in section VI. The 
last section gives a summary and outlook. 
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FIG. 1: A multi-site system is studied in which the two 
outermost sites are attached to electronic reservoirs in ther- 
mal equilibrium with their respective Fermi energies Ef,i and 
Ef, t - The coupling of the outermost sites to the correspond- 
ing lead is described by a spectral density function Jr(lj) and 
the different sites are connected to each other by a hopping 
element A. The on-site energies Ei of the wire can be ma- 
nipulated with a time-dependent external electric field U(t). 



II. MODEL 

The system of interest can be represented in a very 
general form by a time-dependent Hamiltonian 

H(t) = H s (t)+H R + H S R, (1) 

which consists of the time-dependent part describing the 
relevant system Hs(t), a Hamiltonian describing the en- 
vironment Hr, in the presented study the electronic 
leads, and a coupling term H$r between the relevant 
system and its environment. The wire is formed by elec- 
tronic sites n coupled to each other by a hopping element 
A. In second quantization, this orbital tight-binding de- 
scription of the molecular wire reads 

H S (t)=J2 H nn>{t)4 l C n , (2) 
nn f 

where c„ annihilates and creates an electron at site 
n with the anticommutator [cjj,cv] + = 8 n ,n' ■ Here, the 
time-dependence of Hs(t) occurs only in to the irradia- 
tion of the system by an external electromagnetic field 
that manipulates the on-site energies E n with a time- 
dependent on-site potential U n (t). By neglecting possi- 
ble influences of the external field on the tight-binding 
hopping element A, the tight-binding matrix elements 
H nn i (t) can be decomposed into 

H nn ,(t) = -A(8 n+hn/ +S n ^ +1 ) + {E n + U n (t))S nn ,. (3) 

The environment of the wire consists of two electronic 
leads that are modeled by two independent electron reser- 
voirs in thermal equilibrium. For each lead, the Hamil- 
tonian Hr in second quantization is given by 

H R = Y^ u q c\c q (4) 

9 



with c q and c q creating and annihilating an electron in the 
corresponding reservoir mode \q) with mode energy w q . 
Due to the assumed thermal equilibrium of the electronic 
leads, the occupation expectation values of the reservoir 
modes are determined by 

( c l c q') = n F(u q ~ E F )5gqi, (5) 

where rip is the Fermi function and Ep the Fermi en- 
ergy. The left electronic lead is coupled to the first site 
and the right one to the last site N of the wire. To keep 
the notation simple, we will only refer to the left lead in 
further derivations but the formalism has to be applied 
to the right lead as well. The coupling of the left elec- 
tronic lead with the first site of the wire reads in second 
quantization 

H S R = Y J {V q c\c q + V;c) q c 1 ) (6) 

with a system-lead coupling strength V q for each reser- 
voir mode. In general, these coupling values are deter- 
mined by the electronic bands of the leads, e.g. the gold 
contacts, that couple with the energy levels of the wire 
and by the occupation level of these bands with electrons 
given by the Fermi energy and the corresponding Fermi 
function. 



III. REDUCED DENSITY MATRIX AND 
TIME-NONLOCAL APPROACH 

A. Equation of motion 

In general, the equation of motion (EOM) for the com- 
plete density operator p(t) including the wire and the 
leads is given by the Liouville equation 

^ = -l[H{t),p{t)] = - l -C{t)p{t). (7) 

But the information of interest is limited only to the sys- 
tem part of the density operator psif), i.e. the wire part, 
which can be obtained by defining a projection opera- 
tor P projecting the complete system onto the relevant 
part and by tracing out the reservoir degrees of free- 
dom Pp(t) — pRtiR{p(t)} — pr® ps{t)- The irrelevant 
part can be determined by the projection Q = 1 — P. 
Solving the EOM for the irrelevant part to express the 
EOM of the relevant part leads directly to the Nakajima- 
Zwanzig^&SLii operator identity 

Pp = -iPC{t)Pp{t) (8) 
-PC(t) f At l fe-^'' ATQC{T) QC{t')Pp{t') 1 

where the initial-value term is neglected. Since all the 
operators are chronologically ordered in time, the liter- 
ature refers to this TNL approach often as chronologi- 
cal time ordering prescription (COP) 42,43,44 or time con- 
volution approach 4 ^. In order to derive an applicable 
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method to solve Eq. J3J), we apply second order pertur- 
bation theorjs&iSiH and trace over the reservoir degrees 
of freedom to get 

p s {t) = -iC s {t)ps{i) 

-tr R ^£ SR {t)£ dt'U (t,t')£ SR (t')p(t')} (9) 

Here £$, £r and £sr. are the Liouville operators which, 
respectively, apply Hs, Hr and Usr- Uo(t,t') — 
Te~% Jf dT l^s(r)+c R ] j g -f^g time evolution operator of the 
full system without the interaction part H$r and T is the 
time ordering operator in positive time direction. The an- 
nihilation and creation operators of the system and the 
environment are defined in different Hilbert spaces which 
makes it possible to rewrite the coupling Hamiltonian in 
Eq. iJBJ as 

2 

H SR = Y,K X ®$ X (10) 

x=l 

with = E q V q c q , $2 - T, q V*4 and K x = c\, 
K2 = c\. Inserting Eq. IjlOl) in Eq. and using the 
definition of the Liouville operators, a non-Mar kovian 
quantum master equation in the Schrodinger picture is 
obtained 19 . 20 ! 39 . 48 

Ps (t) = -i£ s (t)ps(t) 

"* *<> XX' 

x[K x ,$ x ,,p s (lf,to)]e iH «l t -V]} , (11) 

with U s (t,t') = T + e-iti dTC s(r) being the time evo i u _ 
tion operator of the relevant system. 

The trace over the lead degrees of freedom and all the 
reservoir operators in the dissipation term can be sum- 
marized into correlation functions 

Cxx> (t) = tr R {e iHRt $ x e- iHRt § x ,p R } (12) 

that contain all the information about the reservoir and 
its interaction with the corresponding wire site. For the 
system of interest, these correlation functions decay in 
time thereby causing a memory loss in the dissipation 
term of Eq. ifTTjl. 

Due to the thermal equilibrium condition of the elec- 
tronic leads, two of the four functions are zero, Cn = 
C22 = 0. This reduces the summation over x and x' to 
the pairs (xx 1 ) = (12) and (xx') = (21). Using the prop- 
erty of the correlation functions that C xx <(t) — C* x ,(—t), 
we can define the following auxiliary operators 

A xx ,{t)= f dt'C xx ,{t-t')U s {t,t')K x ,p s {t'), (13) 

Jt 

K x ,(t)= [ dt'C* x , x (t-t')U s (t,t')p s (t')K x , (14) 

Jt 



to simplify Eq. Ijlip. These auxiliary operators incor- 
porate the memory of the system and weight the time- 
dependent electron transfer between the wire and the 
lead. Expressing the quantum master equation (|11|) in 
terms of the auxiliary operators, we get the final expres- 
sion for the master equation of the reduced density op- 
erator at time t 

p s (t) = -i£s(t)p s (t)-^ K ^ & xx ,(t)-A xx ,(t)]. (15) 

xx' 

B. Spectral decomposition 

In analogy to methods recently developed for systems 
coupled to a bosonic bathi 9 ^ 2 ^ 3 ^ 3 ^, we develop equations 
of motions for the auxiliary operators defined by Eqs. (|13l 
I14|l using a numerical decomposition of the spectral den- 
sity Jr(lo) to decompose the reservoir correlation func- 
tions. Solving the trace in Eq. and making use of 
Eq. JSJ, the nonvanishing correlation functions read 

C12W - \ V i\ + E F )e~^ t (16) 

Q 

C 2 i(t) = ]T \V q \ 2 n F (uj q - E F y^ 1 . (17) 
q 

All the external properties of the fermionic reservoir are 
described by a single quantity, namely the spectral den- 
sity Jr(lu), which can be generated by a superposition of 
weighted S functions 

Jr(lu) =Y,K\V q \ 2 S(io-u; q ). (18) 

Equation (|18|) becomes a smooth function for a dense 
spectrum of the reservoir modes. Using the properties of 
the S function to transform the summation over q into 
an integral over lj and extending its lower limit to —00 
by assuming Jr(uj) = for u> < 0, since we need an im- 
proper form of the integrals to solve them by applying the 
theorem of residues, we get the final integral equations 
for the reservoir correlation functions 

f°° du 

C 12 {t)= —j R (w)n F {-u q + E F )e- %UJt . (19) 

J-00 7T 

/OO J 
— J R (w)n F (cu q ~E F )e^ t (20) 
-00 7T 

To solve these integrals, we will pursue a method pro- 
posed by Meier and Tannor— for bosonic systems and 
extend it to fermionic systems using a numerical decom- 
position of the spectral density for u> > 

m 1 

J ^) = EjFT7 o ^r2 > ( 21 ) 

40 fc [w - S2 fe ) + Tf. 
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with real fitting parameters pk, fife and Tfe. This decom- 
position is not restricted to a certain shape of the spectral 
density and can therefore be used to approximate com- 
plicated band structures. This enables one to avoid the 
assumption of the wide-band limit and to take influences 
of the band structure on the dissipative electron transfer 
between the wire and the leads fully into account. 

With the complex roots of the Fermi function and of 
function J2J, the theorem of residues applied to Eqs. 
(IH 120)) results in 

m 
k—1 

9? m ' 

-jYsM^e'^ (22) 

m 

fc=i K K 
2? m ' 

-^E J «K) e " fct ( 23 ) 

with the abbreviations Clt = f2fc+ir/. and fl^ = fife — iFfe 
and the Matsubara frequencies i/fc given by Vk — i 27r ^ +7r + 
These equations determine the coefficients necessary 
to finally write the correlation functions as a superpostion 
of weighted exponential functions: 

m+m' 

Cn{t)= £ a* 2 e^ 4 (24) 

fe=l 

m+m' 

CaiW = E ^i^ 1 '- (25) 
fc=l 

Rigorously, the sum over the Matsubara frequencies 
would be infinite but it can be truncated at a finite value 
depending on the temperature of the system T and the 
spectral width of Jr(oj). This representation for the cor- 
relation function allows us to derive a set of differential 
equations for the auxiliary density operators (|13II14[) . viz. 

^A^(t) = a k xx ,K x ,ps(t)-i[H s (t),A k xx ,(t)] 

+7L'AL'W, (26) 

^A^(i) = (a k , x )* p s (t)K x , -i[H s (t),l k xx ,(t)} 

+ (7* fe 'x)*A^(*) (27) 

with A-xx' (t) = Er=r' A xx'(*) and A«/(t) = 

EkS™ A^ft). These equations can be solved numer- 
ically using a simple Runge-Kutta method without the 



need of diagonalizing the Hamiltonian and together with 
Eq. (TSJ, one now has a complete set of differential equa- 
tions describing the population dynamics in the wire cou- 
pled to the lead in second-order perturbation theory for 
an arbitrary time-dependent wire Hamiltonian. Regard- 
ing the right lead, one just needs to add a second dissi- 
pation term to the master equation l|15|l with differently 
defined K x operators, i.e. acting on the last wire site N, 
and a second corresponding set of differential equations 
for the auxiliary operators. 



IV. CURRENT EQUATION 

An intuitive approach to the electric current equa- 
tion is to consider the rate of change of the num- 
ber of electrons, with elementary charge e, in- 
side the lead which reads in the density matrix 
fjrmalic m 17 - 18 - 19 - 20 - 24 - 26 - 27 - 49 - 50 - 51 - 52 - 53 

h(t) = e^tv{N lP (t)} = -ietv{[N h H(t)]p(t)}. (28) 

Here Ni = c q c q denotes the electron number opera- 
tor of the left lead with the summation performed over 
the reservoir degrees of freedom. Similar to the last sec- 
tions, all calculations refer only to the left lead but are 
also valid for the right lead by adding the corresponding 
terms to the final set of differential equations. The trace 
and the density operator p in Eq. 128|) arc defined in the 
Hilbert space of the full system consisting of wire and 
lead. Solving the commutator in Eq. (|28|l and making 
use of the equilibrium condition of the fermionic reser- 
voir results in Ii(t) = — 2e l m ti{V q c\cqp(t)}. Since 
the trace is invariant under a transformation into the in- 
teraction picture, A(t) = Us{t,to)A for any operator A 
defined in the Hilbert space of the full system, one can 
use the integrated form of the Fiouville equation in the 
interaction picture p(t) = p(to) — i L dt' [HsR(t') t p(t')] 
for the time evolution of the full density operator p(t) 
in the interaction representation. After solving all the 
commutator relations, separating the system part from 
the reservoir part and employing again the definition of 
the correlation functions (|12|) . the ansatz for the current 
iPH finally reads 

J,(t) = 2eRe f dt' (tr s {ct(i)ci(i')ps(*')} C 12 (t - t') 

-tr s {ci(f)^(*)ps(0} C3i(*-*0) • (29) 

The current information is partly contained in the tempo- 
ral phase relation between the annihilation and creation 
operators which can be incorporated into a single time 
evolution operator Us, defined in the Hilbert space of 
the relavant system, by using the cyclic properties of the 
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trace to write Eq. I|29|) as 
h{t) = 2eRe (tr g < r 
-tr s 



with the modified corresponding auxiliary operators 



l , dt'UHt'^^psit')) C 12 (t-t') 

to 



dt/C^t-^UsMK* 



(33) 



to 



<4 jT dt'C/t(t',t) Ps (i')ci| Cl^t-t')) • (30) 



These integrals have the same structure as the auxiliary 
density matrices defined in the last section and can be 
expressed in terms of A\2(t) and Ai2(i). Thus, we get 
the final equation for the time-dependent current 

/,(<) = 2eRe (tr s {4A 12 (i) - c\A 12 (t)}) (31) 

between the left lead and the first site of the wire. Similar 
to the master equation l|15l) for the population dynam- 
ics, all the information about the interaction of the sys- 
tem with the reservoir is contained in the time-dependent 
auxiliary matrices (| 1 31 114f> weighting the corresponding 
system operator. Eq. H31(l is valid within the framework 
of assumptions we have made to derive Eq. (|15|) and al- 
lows one to measure the time-resolved current for arbi- 
trary time-dependent systems. 

The EOM for the required auxiliary operators are al- 
ready known from the time propagation of the popula- 
tion. Thus, we have a complete set of differential equa- 
tions consisting of Eqs. I|2t)ll27|) . the master equation 115|) 
and the current equation (|31|l and therefore an applica- 
ble formalism to compute the current and the popula- 
tion dynamics for a wide range of systems. The same 
set of equations has to be applied to the last site of the 
wire as well to describe a wire enclosed by two leads, as 
shown in Fig. ^ Furthermore, the ability to define any 
desired number of leads (with respect to the given com- 
putational resources) makes it possible to realize driven 
current switches with an arbitrarily complicated config- 
uration. 



V. ALTERNATIVE TIME-LOCAL APPROACH 

Since we started with the Nakajima-Zwanzig iden- 
tity to derive our equations for the population dynam- 
ics the final expressions are in a TNL form. In this 
section, we shortly present their TL counterparts with- 
out giving an extended comparison of both approaches. 
In the literature, the TL approach is also known as 
time-convolutionless formalisms 5 -, partial time ordering 
prescriptioni24iMi or Tokuyama-Mori approach 54 . Us- 
ing the Tokuyama-Mori identity one gets in second-order 
perturbation theory for the wire-lead coupling the follow- 
ing equations for the population dynamics of the wire, 



p s {t) = -iC s {t)p s {t) 



(32) 



A xx ,(t)= dt'C*, x (t-t')Us(t,t')K x ,. (34) 

Jt 



Ps(t)A X x'(t)} 



One can also derive these equations by just applying the 
substitution p{t') = Ug(t, t')p(t) to the corresponding set 
of TNL equations^*, neglecting the influence of dissipa- 
tion during the time propagation of the density operator 
within the integral. Thus, the current equation readsi2i2£ 



h{t) = 2eRe (tr s {c\k 12 {t) Ps {t) - c\p s (t)A 12 (t)}^ . 

(35) 

TNL and TL approaches for similar systems have been 
compared in recent numerical studies 39,55 and it depends 
strongly on the coupling parameters which of both is fa- 
vorable. Within the parameter regime we used for the 
following numerical investigations, TL and TNL calcula- 
tions showed almost identical results. 



VI. NUMERICAL RESULTS FOR THE LASER 
DRIVEN MOLECULAR WIRE 

The system depicted in Fig.^represents a simple con- 
figuration which allows one to study a variety of trans- 
port phenomena. The following numerical results are ex- 
pressed in terms of the tight-binding hopping parameter 
A, where A = O.leV is a reasonable value for molecular 
system a 17 ! 24 ! 26 ! 27 ! 50 ! 51 ! 52 ! 53 . The system temperature is 
set to T = 0.25A/fc B = 290K. Despite the fact that it 
would be possible in Eq. (|21|l . we do not simulate a re- 
alistic spectrum for the coupling between the leads and 
the wire, but rather restrict the spectrum to a single 
Lorentzian. For example, a realistic gold s-band would 
be a superposition of different Lorentzians forming cer- 
tain band edges. Taking only one Lorentzian into account 
corresponds to a rough approximation of only one band 
edge, but this makes the system and the underlying phys- 
ical processes easier to understand. Due to the weak cou- 
pling requirement, the spectral density should be about 
one order smaller than the internal dynamics, thus the 
peak of the Lorentzian Jr(lo) is set to 0.1A, what is guar- 
anteed by the condition pi = 0.1 (4AriiT 2 ). Reasonable 
values for the bandwith parameters Tk are in the region 
of lOeV. With the chosen energy settings, a time unit in 
the system corresponds to 0.66fs what enables one to re- 
solve time-dependent effects on a femtosecond scale. The 
resulting current unit can be extracted from Eq. i|31|) and 
is equal to a macroscopic value of l[e] = 2.43 * 10~ 4 A. 

As a basis set for the tight-binding Hamiltonian we 
choose the many-body Fock space in which each site is 
represented by a two-state vector \xn)- The set of these 
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state vectors forms the total state vector by the ten- 
sor product \ty) = xi) ® IX2) ® ■■■ <8> \xn)- The corre- 
sponding annihilation and creation operators for site n 



c f = 





1 



1 




(36) 



represent a complete operator basis for the master and 
the current equation. 



A. The transient oscillation of the undriven wire 

The equilibrium condition of the on-site population of 
the wire in the absence of a second lead is determined by 
the Fermi function n F taken at the on-site energy of the 
site, e.g. tr(c\cips(t)) — n F (Ei — Ep). In this particular 
case, the current drops to zero after the wire is filled to its 
equilibrium value, simply due to the lack of a closed elec- 
trical circuit. The situation is more interesting in the case 
where two leads are coupled to the wire, as it is shown in 
Fig-ffl Fig. [3 shows a situation for a wire with two sites, 
in which the wire states at time zero are unoccupied. 
The equilibrium state is reached after a relaxation time 
that mostly depends on the wire-lead coupling strength. 
While the number of electrons in the relevant system is 
not constant in time, the trace of the reduced density op- 
erator is conserved and normalized. A bias voltage on the 
system can be simulated by setting a difference between 
the left and right Fermi energy, here E F ^ — Ep^ = 10 A. 
Since we set the on-site energies to E\ = E2 and apply 
a weak coupling scheme where the internal dynamics of 
the inner system is about one order faster than the elec- 
tron transfer between the leads and the wire, a spatial 
drop of the population within the wire in its equilibrium 
state cannot be observed in the upper panel of Fig. [2] 
The equilibrium population of the entire wire is deter- 
mined by the average over the expectation values of the 
left and right sites, which is given by tr(c\aps(t)) — 

tr(4caps(t)) = {n F (E x - E Fl ) + n F {E 2 - E Fr ))/2 if 
the spectral densities of both leads, i.e. their coupling 
to the wire, are the same. The transient oscillations of 
the corresponding directed net current decay and a con- 
stant equilibrium value is reached, when the relaxation 
process is over. The upper bound for the net current 
is determined by the small wire-lead coupling strength, 
namely the spectral function Jr(uj) with its maximum 
value of 0.1A. 



B. AC voltage 

Time-dependent electronic reservoirs in molecular sys- 
tems can be modelled within the Floquet theory by ap- 
plying a gauge transformation on the wire Hamiltonian 56 
or by using a source-Redfieldi^ equation. Despite the 
fact that we initially assumed that the lead Hamiltonian 
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FIG. 2: Upper panel: population dynamics measured by 
tr (CnC„ps(t)) (n — 1,2) for the system with two sites with- 
out driving field U(t) as a function of time t starting with 
unoccupied wire states. The applied DC voltage is realized 
by a difference between the left and the right Fermi energy 
of Ep t i — Ef, t = 10A. The on-site energies of the wire are 
aligned and centered between the left and the right Fermi en- 
ergy E\ = E2 — Ef,t + 5A = E F .i — 5A. Lower panel: the 
corresponding currents flowing from the left lead into the first 
site heft, from the right lead into right site I right and the net 
current given by I net = {heft - Iri g ht)/2. 



is time- independent, the formalism allows one to treat 
time-dependent lead Hamiltonians as well by simply con- 
sidering time-dependent coefficients for the correlation 
functions C(t) = J2k a k(t)e lk ^ t , in the case of the TL 
theory, and a time derivative given by 

, m+m' 



-c(t) = MtW k{t)t hk(t) + j k (t)t) 

k 

m+rn 

+ «fcWe 7fc(t)t . (37) 



To get the required form C(t) = jk(t)C(t) which is nec- 
essary to derive differential equations for Eqs. 1131 114fl . 
we have to assume that 



\a k {t)\ « |ofe(*)7fc(*)l, 



? lk{t)t jk(t)t\ « |e 7fc(t) *7fcWI 



(38) 



(39) 



for all times t not negligible in the time integration of 
the auxiliary operators in Eqs. (|13l I14fl . In general, 
this assumption corresponds to the statement that the 
lead dynamics is sufficiently slower than the lead-wire 
coupling dynamics. The lead- wire coupling can be in- 
fluenced by a variation of the chemical potential of the 
electron distribution in the leads, by changing the param- 
eters fifc, Tk,Pk of the spectral decomposition, or both. A 
change in the chemical potential and keeping the spectral 
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FIG. 3: The AC voltage (upper panel) with amplitude Vb and 
frequency uiv = 0.00625A drives an AC I[e] (bottom panel) 
between the left and the right lead through an undriven and 
not disordered wire consisting of five sites. A current cut-off 
appears for large AC voltages, since the wire-lead coupling is 
weak. 



function unchanged would correspond to a charging pro- 
cess, while shifting the spectral density, i.e. increasing or 
decreasing the Qk parameters, and the chemical potential 
by the same value describes a total change of the electro- 
static potential without an alteration of the population 
in the leads. In the present work, we restrict ourselves to 
a time-periodic modulation of the difference between the 
chemical potentials of the left and right lead Ep^—Epj = 
Vq sin(wyi), thereby approximating an AC voltage on the 
leads with amplitude Vb- In this special case, the quality 
of the approximation of slow dynamics can be determined 
considering the general time dependence of the coeffi- 
cients given by Eqs. <E21 E3J as |a fe (t)| ~ \n F (E F (t)) \ 
for the non-Matsubara terms with the corresponding 
time-independent exponents "fk- Therefore, the time 
derivative becomes \ak{t)\ ~ \n% (E F (t)) (3E F (i)\ ~ 
l^^cos (uj v t) n 2 p (E F (t)) | and conditions J3H1 EU are 
justified if Va £ v | <C 1. Regarding the Matsubara terms, 
the relevant time-dependent coefficients are jk(t) = 
-(2nk + ?r)T + iE F (t) and \a k (t)\ = \2TJ(v k )\ for k = 
1 ... ml and it can be shown that conditions (pIHl are 
also fulfilled if (^zntyT ^ ^ which complies with the 
condition for the non-Matsubara coefficients. A disregard 
of the conditions (|38l I39|) would cause additional fast os- 
cillations of the current during those times in which the 
AC voltage changes rapidly. The AC voltage causes an 
AC in the system, both shown in Fig. El that follows the 
driving voltage. For large voltage amplitudes Vb a finite 
cut-off appears due to a bottle neck of the system given 
by its weak coupling to the electronic leads. 



C. Coherent destruction of tunneling 

Coherent destruction of tunneling 24 i 26 i 52 i 53 i 57 (CDT) is 
a well understood quantum mechanical effect where an 
external periodic driving field of the form 

A(t) = A sin (ujt) (40) 

yielding an asymmetric potential U n (t) = A{t)5\ n — 
A(t)S2n in Eq. @ can suppress the time-averaged cur- 
rent in the driven two-site system depicted in Fig. ^ 
The current completely breaks down if the ratio of the 
field parameters Ao/fuvd is equal to a zero of the Bessel 
function J (e.g. 2.405, 5.520, 8.654, ...), a condition 
that holds for isolated quantum systems as well as for 
open systemsS^S^SiS. This is a simple but nontriv- 
ial system for a time-dependent conduction formalism 
and we extend the results of former calculations which 
applied the wide-band approximation^SSiSiS to finite 
band effects by using a spectral decomposition. The fi- 
nite width of the decomposed spectral density (|21|l causes 
an additional decay of the current with increasing ampli- 
tude since the coupling strength between the lead and 
the corresponding site decreases when the on-site en- 
ergy of the coupled site is driven to the edges of the 
Lorentzian during the oscillation generated by the exter- 
nal field A(t). Naturally, this effect becomes more dom- 
inant for a smaller bandwidth parameter Ti, shown in 
Fig. The time-averaged current shows the predicted 
breakdowns at the zeros of Jo- The amplitude Aq of the 
applied electrical field is measured in dimensions of [A] 
which corresponds to field strengths of about 10 8 ^ by 
assuming intra-atomar distances of lA. The driving fre- 
quency is set to LOd — 10A which is in the low energy 
branch of infrared light. Here the bias voltage is set to 
E F i - Ep. r = 60A = 6 eV. 

D. Optical current switching 

The major advantage of the derived conduction formal- 
ism, in addition to the avoidance of the wide-band limit, 
is the applicability to unrestricted time-dependent sys- 
tems. Former CDT studies of currents in open quantum 
systems^* 2 ^^ were based on an infinite-time averaging 
of the currents due to the mathematical nature of the 
used approaches. In Fig. we apply a finite laser pulse 
with a Gaussian shaped amplitude 

A{t) = Ao exp ) sin (cu d t) (41) 

to the asymmetric driven system described in the last 
subsection. The peak amplitude of the Gaussian was set 
to A a = 24.05A, where the CDT relation applies. This 
finite laser pulse causes the time-resolved current to os- 
cillate around zero, shown in the center panel of Fig. [S] 
The situation becomes more obvious by looking at the 
time averaged current depicted in the bottom panel of 
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A [A] 

FIG. 4: The time-averaged current for the periodically driven 
two-site system as a function of the amplitude for different 
bandwith parameters Ti. The frequency of the external field 
is ujd = 10A. The applied bias voltage is 60A. The Fermi 
energies Ef, t = Oi — 30A, Ep,i = fii + 30A and the on-site 
energies E\ = E2 = fii are set in relation to the parameter 
fii of the spectral density. 



30 -i 




FIG. 5: The upper panel shows the Gaussian excitation pulse 
with a peak amplitude of Ao = 24.05A and a width of a 2 — 
50. It excites the time-resolved net current depicted in the 
middle panel. If the net current is averaged over three periods 
of the fast pulse oscillation, a complete suppression can be 
observed, depicted in the bottom panel. The on-site energies 
are equal E\ — E2 = fii and centered between the left and the 
right Fermi energy, which define a DC voltage of Ef, t — Ep,i = 
I0A. 

Fig. where the averaging was taken over three oscilla- 
tions of the driving field. The averaged current is almost 
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Phys. 93, 77 (1996). 



suppressed during the time in which the Gaussian laser 
field is applied to the system and CDT is performed with 
a finite laser pulse. The effect is superimposed by tran- 
sient oscillations because the system is under the con- 
straint of permanently changing variables and tries to 
find its equilibrium state. This finite-time effect gives 
rise to interesting experimental realizations and further 
investigations will deal with the idea to apply the cur- 
rent formalism to non-local optimal control algorithms 
to compute more sophisticated and applicable external 
control fields. 



VII. SUMMARY AND CONCLUSION 

We have developed a time-dependent non-Markovian 
conduction formalism based on a projection operator 
approach for the density matrix with a numerical de- 
composition of the spectral density which allows one to 
study time-resolved currents and population dynamics in 
molecular wires for arbitrary time-dependent wire Hamil- 
tonians, which is an useful extension to existing theories 
and applicable for a wide range of systems. The for- 
malism includes the coupling to the electronic leads in 
second-order perturbation theory and a non-perturbativc 
interaction with external fields. Its validity was shown 
for the example of coherent destruction of tunneling in 
a driven two-state system and further effects like optical 
control of current by using a short laser pulse, alternat- 
ing currents and electrical relaxation of the system into 
a biased equilibrium state were presented. These effects 
were investigated on a femtosecond time scale, which 
is an important aspect for the the applicability of the 
formalism for future investigations of femtosecond spec- 
troscopy, optical control and optimization of electrical 
currents. One can easily include electron- vibrational cou- 
plings in second-order perturbation theory and electron- 
electron correlations resulting in a powerful tool to sim- 
ulate currents under quite realistic conditions as long as 
the weak coupling scheme is justified. Furthermore, re- 
cent publicationsi2i2£ suggest a way to avoid this weak 
coupling condition within the framework of the presented 
formalism. 
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